Numerical study of the critical behavior of the Ashkin- Teller model at a line defect 
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We consider the Ashkin- Teller model on the square lattice, which is represented by two Ising 
models (cr and r) having a four-spin coupling of strength, e, between them. We introduce an 
asymmetric defect line in the system along which the couplings in the a Ising model are modified. 
In the Hamiltonian version of the model we study the scaling behavior of the critical magnetization 
at the defect, both for a and for r spins by density matrix renormalization. For e > we observe 
identical scaling for a and r spins, whereas for e < one model becomes locally ordered and the 
other locally disordered. This is different of the critical behavior of the uncoupled model (e = 0) 
and is in contradiction with the results of recent field-theoretical calculations. 



I. INTRODUCTION 

In a system which is divided into two parts by a defect 
plane translational invariance is broken and the physi- 
cal properties are different in the defect region, which 
has a width of the correlation length, ^. At the critical 
point where ^ is divergent the scaling properties of local 
quantities, such as the defect magnetization or the spin- 
spin correlation function, could be different from that in 
the bulk[l]. Relevance or irrelevance of the perturbation 
caused by a weak defect can be analyzed within the frame 
of phenomenological scaling theory [^-|5|- the perturba- 
tion is coupled to the local energy operator, then the bulk 
fixed point is stable if the correlation length critical ex- 
ponent of the pure system, v, is greater than 1. In this 
case the defect exponents are the same as at any other 
point of the bulk. In the opposite case for v < I, gener- 
ally a new fixed point governs the local critical behavior, 
the properties of which are different for weakened and for 
enhanced local couplings. For weakened defect couplings 
the defect usually renormalizes to a cut and the local 
critical behavior is the same as at the ordinary surface 
critical point On the contrary for enhanced defect 

couplings the defect usually renormalizes to an ordered 
interface and the local critical exponents are the same as 
at the extraordinary surface transition [6H8|. Examples 
for modified defect critical behavior can be found in the 
two-dimensional (2D) q = 3 state Potts modelj^, [l^], in 
the Baxter- Wu model [H] or in the Ashkin- Teller (AT) 
modelH. 

According to scaling considerations in a system having 
a correlation length exponent v — I a defect is a marginal 
perturbation. This happens for the 2D Ising model for 
which the local magnetizatio n. 11 3l and the spin-spin cor- 
relation function at the defect [Tj] has been exactly calcu- 
lated and the local critical exponent of the magnetization 
is found to be a continuous function of the strength of the 
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defect. These results have been generalized recently for 
inhomogeneous line defects, in which the defect coupling 
is a smooth function of the positionfl5|. 

In a system which has more complicated interaction 
between the spin variables one can define special defects 
which influence only some parts of the interaction. This 
type of asymmetric defect has been considered recently 
by NaonJJj in the ATpT^ and in the Baxter models [l^]. 
We remind that both the AT and the Baxter models are 
expressed in terms of two sets of Ising spin variables [Tj^, 
say a and t, with two-site {K2) and four-site {K4) in- 
teractions. The asymmetric defect introduced by Naon 
is represented by a line of modified two-spin couplings 
in one of the Ising models (say for the a spins), whereas 
the other interactions are left unchanged. Using field- 
theoretical methods Naon has determined the decay of 
the critical spin-spin correlation function, both for the 
a and for the t spins and the calculated local magne- 
tization exponents are found to be independent of the 
interaction between the Ising models, which is measured 
by e = Ki/K2. Similar results are obtained in the case of 
an asymmetric and inhomogeneous line defect in the AT 
modelfm. These results are somewhat surprising, since 
in the pure AT and Baxter models e is a marginal pertur- 
bation and several (bulk and surface) critical exponents 
are continuous functions _20] of e. 

In this paper we are going to revisit this type of prob- 
lem and study numerically the local critical behavior at 
an asymmetric defect in the AT model. Our aim is to 
shed some light to the physical mechanism which is be- 
hind this problem and to confront the measured criti- 
cal exponents with the results of field-theoretical calcu- 
lations. Here we shall consider both chain and ladder 
type defectsil| and the system is studied by the trans- 
fer matrix, which is perpendicular to the defect. We use 
the Hamiltonian limit [21', l22'| of the transfer matrix and 
calculate the local magnetization both for a and t spins 
by density matrix renormalization[2^ (DMRG). At the 
critical point the local magnetization exponents are then 
deduced through finite-size scaling. 

The structure of the paper is the following. The model 
and its basic properties are described in Sec. [Hi Numer- 
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ical results of the scaling of the interface magnetization 
are presented in Sec. Illll and discussed in Sec. HV] 



II. THE AT MODEL WITH A DEFECT LINE 

The AT model is defined in terms of a four-state spin 
variable [TtI. which is represented by a pair of Ising spins: 
CTi = ±1 and Ti = ±1 at the lattice site i. Between the 
same set of spins there is nearest neighbor interaction 
of strength K2 and the two Ising models are coupled by 
the product of the energy densities, which is given by a 
four-spin term: Kj^UiGjTiTj, where i and j are nearest 
neighbors. We consider the system on the square lattice 
and work with the row-to- row transfer matrix Tat- In 
the Hamiltonian limit the transfer matrix can be written 
as Tat ^ exp(— kHat), where k is the lattice spacing 
in the "time" direction and "Hat is a one-dimensional 
quantum Hamiltonian given bvj2li [2^ 

■hat= - E«'^m+<^i)-'^E«+^) 



Here cr^'^ and t^'^' are two sets of Pauli matrices at site i 
and h is the strength of the transverse field, which plays 
the role of the temperature in the classical system. As be- 
fore the ratio of the couplings is denoted by e = K^jK^- 
The system in Eq.(IT]) is self-dual and the self-duality 
line: he = I represents the critical line separating the 
ferromagnetic and the paramagnetic phases of the system 
for -1/^2 < e < 1. In the region -1 < e < -1/^2, 
there is a so-called "critical fan" , which extends to both 
sides of the self-duality line and in which the system stays 
critical |21i] . 

The critical properties of the AT model are exactly 
known through conformal invariance and Coulomb-gas 
mapping 24, .25j . The decay of the spin-spin correlations: 
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is characterized by the anomalous dimension, Xm = 1/8, 
which does not depend on the value of the coupling e. On 
the contrary the decay of the connected energy-energy 
correlations: 
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involves the anomalous dimension: Xe — 
7r/[2 arccos(— e)], which is e dependent [2l| . Similar 
holds for the correlation length critical exponent z/, 
which is given by: V = 1/(2 - Xe) for -1/^2 < e < 1 
whereas it is formally infinite in the critical fan. 

Finally, the decay of the end-to-end correlation func- 
tion at the critical point involves the surface magnetiza- 
tion scaling dimension, x'^: 



which is also coupling dependent [26| : x"^ = 
arccos(— e)/7r. 

The marginal operator for the AT model is associ- 
ated with the four-spin term, erf (Tf_|_j^r/r/_^j^, which has 
a scaling dimension, Xi^ = 2 = D, independently of e, for 
-1< e < 1. 



A. Ladder and chain defects 

A line defect in the 2D classical model is put in the 
"time" direction and can be of two types: chain defect 
or ladder defect (see Fig. 6.1 of Ref. [l|). In the prob- 
lem studied by Naon 16] and we consider here a line of 
two-spin couplings between the cr-spins are modified. In 
the Hamiltonian limit, when the defect is placed between 
sites i = and i — \ (ladder defect) the perturbation is 
given by: 
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where J is the strength of the defect. On the other hand 
for a chain defect with modified two-spin couplings at line 
« = the perturbation in the Hamiltonian limit involves 
the term: 
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The ladder and chain defects transforms into each other 
through duality [l3| and their strengths are related as: 
li o 1/J. 

For two decoupled Ising models with e = the local 
magnetization exponents at the defect site, i = 0, gener- 
ally are different for the cr and for the r spins, which are 
denoted by a;^ and a;J^, respectively. While x^^ keeps its 
bulk value, is a continuous function of the strength 
of the defect [li, Hi]: 
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for a ladder defect and 
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0, (7) 



6 = 0, (8) 
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for a chain defect. For a spins the marginal operator is 
the local energy density, which has its anomalous dimen- 
sion = 1, independently of the value of J or h. 

If we switch on the interaction between the Ising mod- 
els the defect critical behavior could be modified. For 
small e the anomalous dimension of the four-spin opera- 
tor is x^ — x1~\-x],— 2, which is just the marginal value. 
We note, however, that field-theoretical calculations of 
the critical spin-spin correlation function [l6{ come to the 
conclusion, that e is an irrelevant variable. Here we re- 
visit this problem by numerical methods. 
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III. NUMERICAL STUDY 

Here we consider the Hamiltonian version of the AT 
model at a finite lattice with —L < i < L and use fixed- 
spin boundary conditions: a^j^ — t^j^ = +1. Using the 
DMRG method we calculate the ground-state expecta- 
tion value of the magnetization operators at the defect, 
(0|crjj|0) = m^{L) and (0|ro^|0) = ml{L). The parame- 
ters in Eq.(IT|) are taken at different points of the critical 
line: h = he — 1 and — 1 < e < 1 and we have considered 
different strengths of the defect: h > and J > 0. Ac- 
cording to finite-size scaling theory these magnetizations 
asymptotically behave as: 



(9) 



In the numerical calculation we went up to L = 85 and 
in the DMRG method we have generally kept around 
m = 150 states of the density matrix in order to obtain 
a good numerical accuracy. 

From the values of the matrix element in Eq. © at two 
different sizes, L and bL, we have calculated effective, 
size-dependent exponents through two-point fits: 



lnm^(&i) -lnmfJ(L) 



(10) 



In order to obtain the same numerical accuracy for the 
different lengths, we keep the ratio b between neighboring 
sizes approximately constant. The effective exponents 
evolve towards their exact values when the mean size 
associated with the two-point fit, (L) — L(b+l)/2, tends 
to infinity. 

In the actual calculation we consider first in Sec lIII Al 
the known limiting cases (decoupling limit, system with- 
out defect) in order to check the accuracy of the numeri- 
cal method. Afterward we study the general model with 
defect for e > and e < in Sec lIII Bl and in Sec lIII C[ 
respectively. 

The numerical calculations are performed on several 
PC-s, the total CPU time being equivalent to ^ 250 days 
in a single 2.4 GHz processor. We note that an alterna- 
tive numerical method is to use Monte Carlo (MC) sim- 
ulations on the classical model with L x M sites. As can 
be seen in an investigation of a related problem|27| the 
DMRG and the MC methods generally provide results 
with approximately the same accuracy using the same 
CPU times. 



A. Decoupled or non-defected systems 

We start with the decoupled model and calculate the 
critical local magnetization of the a Ising model having 
a ladder defect. The finite-size magnetization exponents, 
x'^{L), which have been obtained through Eq. pO)) are 
plotted as a function of 1/L in the left panel of FiglT] 
for different strength of the defect, J. Extrapolating for 
large L the limiting exponents are in good agreement 
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FIG. 1: Numerical estimates of the magnetization exponent 
(L) for different finite systems of length L at special points 
of the phase diagram. Left panel: at the decoupling point, 
e = 0, for various values of the strength of the ladder defect, 
J. Right panel: non-defected AT model for different values of 
the coupling, e. The exact results are indicated by horizontal 
broken lines. 



with the analytical results in Eq.([7]), which are indicated 
by dashed lines. 

We have also studied the AT model without defects, 
but for different values of the coupling, e > 0. The finite- 
size magnetization exponents, which in this case corre- 
spond to the bulk exponents, are plotted as a function 
of 1/L in the right panel of Fig[T] As seen in the figure 
the extrapolated exponents are fairly close to the exact 
value, Xm = 1/8, which indeed does not depend on e. The 
correction terms, however, are larger for a more strongly 
coupled system. 



B. Positive coupling: e > 

This part of the phase diagram of the Hamiltonian 
model corresponds to the classical model with positive 
Boltzmann-weights, thus the results obtained in this do- 
main should be compared with the field-theoretical cal- 
culations in Ref . [l6| . First we present in Fig 12] the finite- 
size magnetization exponents (both for a and t spins) 
which are calculated at the largest length, L — 85. 
Here we have a system with a chain defect of strength, 
h = 1/2, 2/3, 3/4, 1., 4/3, 3/2 and 2 and various posi- 
tive values of the coupling. For the sake of comparison we 
also present the analytical results in the decoupling limit, 
6 = 0, which are plotted with dotted-dashed and broken 
lines for a and r spins, respectively. We note that ac- 
cording to field-theoretical investigations [ll] these results 
should hold for e > 0, too. The numerical results in Figl2] 
seem to be in contradiction with the field-theoretical con- 
jectures in two respects, i) The effective (L-dependent) 
magnetization exponents vary with the coupling, e > 0, 
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FIG. 2: Finite-size magnetization exponents x'^{L) calculated 
at the largest length L at a chain defect of strength h and for 
positive values of the coupling, e > 0. The open (filled) sym- 
bols are for a — a {a = t) spins. Note, the very similar scaling 
behavior of a and r spins. The exact results for the decoupled 
model with e = are indicated by dotted and broken lines, 
for a and r spins, respectively. These represent the field- 
theoretical result Insets: Size-dependence of the effective 
magnetization exponents for a spins, a) For enhanced defect 
strength, h — 2/3, the finite-size exponents are expected to 
approach x^{L) — 0, which corresponds to an EI fixed point, 
b) For reduced defect strength, h — 3/2, the finite-size ex- 
ponents are expected to approach x"^ — x^^ = xf^, which 
corresponds to an OI fixed point. 



in particular for large e there is a considerable difference 
from the value at e = 0. ii) The magnetization exponents 
at the T spins are different from the predicted bulk value, 
Xm = 1/8, and these are close to that measured values 
at cr spins at the same system. The trend of the effective 
exponents with the size of the system is different for en- 
hanced {J > I or h < 1) and reduced {J < 1 or h > 1) 
defect couplings, respectively. 

For enhanced defect couplings the local magnetization 
exponents are smaller than the values at the decoupling 
limit and the finite-size exponents - for sufficiently large 
sizes - are decreasing with L. This is illustrated in the 
inset a) of Figl2]for h = 2/3. We expect that this decreas- 
ing tendency will continue for larger sizes and the extrap- 
olated exponents will approach = = 0. This type 
of behavior is characteristic for an ordered defect (see 
Eqs.iO and ([8|) for J — > oo and h ^ 0, respectively), 
and we assume that the local transition for enhanced 
couplings is governed by an extraordinary interface (EI) 
fixed point. 

For reduced defect couplings the effective, in- 
dependent exponents are larger than their value at 
the decoupling point and these are continuously increas- 
ing with L. In order to check a possible overshooting 
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FIG. 3: The same as in Fig[5]but for negative coupling, e < 0. 
Note, the complementary scaling behavior of a and t spins: 
Xmih) ~ x^(l//i). Insets: Size-dependence of the effective 
magnetization exponents for a spins, a) For enhanced defect 
strength, h = 2/3, the finite-size exponents are expected to 
approach x'^{L) = 0, which corresponds to an EI fixed point, 
b) For reduced defect strength, h = 3/2, the finite-size ex- 
ponents are expected to approach a non-trivial value: x'^{e), 
which corresponds to a special OI fixed point. 



effect we have made calculations up to i ~ 500, but no 
bending of the curves are found. We have also compared 
the magnetization exponents at the two sublattices, 
i.e. x'^ with x^„. For a finite L, x'^{L) is larger than 
a;^(L), however the limiting value of the two exponents 
are very close to each other, see Fig 15] for h > 1. We 
expect, that for the true exponents we have x'^ — x^. 
A more difficult point is to decide about the eventual h- 
dependence of the exponents. Here one should note that 
the effective exponents show a quite strong h, as well as 
L-dependence. As illustrated in the inset b) of Figl2]the 
extrapolation of the data with L is very difficult, at least 
for not too large values of h. To interpret the data we 
use the hypothesis, that the interface critical behavior 
is controlled by an ordinary interface (01) fixed point. 
In this scenario the true exponents are independent of 
h and their values are the same as at a free surface, as 
given below Eq.Q. Indeed, for larger values of h the 
extrapolated exponents are fairly close to xf^{e). 



C. Negative coupling: e < 

This part of the phase diagram of the Hamiltonian 
model corresponds to the classical model with negative 
Boltzmann- weights. The calculated finite-size magneti- 
zation exponents for the same absolute values of e and for 
the same values of h as in Sec lIIIB] are shown in Fig|31 
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The two figures, Figl3]and Figl2l show several differences. 
The most important difference is that here the limiting 
values of and are different. The effective expo- 
nents for a spins with h are approximately equal to the 
effective exponents for r spins, however at 1/h. As far 
as the defect exponents on the a spins are concerned the 
trend is similar as for e > 0. For enhanced defect cou- 
plings the finite-size exponents seem to approach x"^ = 0, 
thus the a spins are locally ordered and their critical be- 
havior at the defect is controlled by the EI fixed point. 
On the contrary for reduced couplings at the defect the 
finite-size exponents are monotonously increasing with 
the size and they seem to approach a non-trivial value: 
ccj^ = Xdof(e). We expect, that this limiting value de- 
pends only on e, but does not depend on J < 1 or ft, > 1, 
thus it can be identified as a special interface exponent 
in the model with a vanishing bond, J = 0. 



IV. DISCUSSION 

The numerical results presented in the previous Section 
about the interface critical behavior of the AT model with 
a (asymmetric) line defect can be interpreted in terms of 
an RG phase phase diagram, which is shown in the right 
panel of Fig|3] both for the a and the r spins. As a 
comparison in the left panel of FigU we show the RG 
phase diagram for such a (symmetric) defect, which is 
proportional with the local energy density and given by: 



Vcha 



(11) 



for a chain defect and similarly for a ladder defect. In 
this case the a and r spins play equivalent role and the 
interface critical behavior has been studied previously in 
Refs.[TTI. [12} . According to these results, which are in 
agreement with a relevance- irrelevance analysis 0-01, the 
bulk fixed point (B) in the system is unstable for e > 0. 
The numerical results also indicate that the critical be- 
havior at the defect is governed by the 01 and EI fixed 
points, for reduced and enhanced strength of the defect, 
respectively. This type of RG phase diagram is suggested 
to be valid for the asymmetric defect too, as illustrated 
in the right panel of Fig|31 This result indicate, that the 
phase-diagram in the left panel of Fig U for e > is prob- 
ably valid in a coarse-grained description: the direction 
of the fiow at B depends only on the condition, if the 
local energy at the defect is reduced or enhanced with 
respect to the pure system. 

For negative coupling, e < 0, the structure of the 
phase-diagrams are changed, both for symmetric and 
asymmetric defects. For symmetric defects the bulk 
fixed-point becomes the stable one[ll|,[l3|. It is certainly 
understandable, that the EI fixed point is unstable, since 
for e < ordered defects would bring a positive contri- 
bution to the energy. The same reasoning holds for the 



asymmetric defect, too. Thus for e < the a and the r 
spins can not be ordered at the same time. In this case 
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FIG. 4: Schematic RG phase diagram at a defect line of 
strength J ~ in the critical AT model. Left panel: sym- 
metric defect, right panel: asymmetric defect for the a and 
T spins. The RG flow is different for different signs of the 
bulk four-spin coupling e. For e > 0, the flow is towards the 
ordinary interface fixed point (OI) when J < 1 and the ex- 
traordinary interface fixed point (EI) when J > 1, both for 
symmetric and asymmetric defects. When e < 0, the flow is to 
the bulk fixed point for a symmetric defect. For asymmetric 
defect the RG flow is different for a and r spins. 



in the energetically favorable state one set of spins is dis- 
ordered (being in the 01 fixed point) and the other set of 
spins is ordered (EI fixed point). This is the situation we 
have found numerically and which is shown in the phase 
diagram in the right panel of FigU) 

The field-theoretical results about the spin-spin corre- 
lation function at the defect line suggest a different RG 
phase diagram[l^. According to this theory, which ex- 
pected to hold for positive Boltzmann-weights (e > 0) 
the coupling term is an irrelevant perturbation and the 
critical behavior at the defect is the same as for the de- 
coupled systems. Our numerical results are in contra- 
diction with this scenario. In this respect we mention 
that in the asymmetric defect problem two potentially 
marginal operators are involved. One is the product of 
the energy-densities in the AT model and the second is 
the local energy-density (line defect) in the Ising model 
with a spins. The interplay of these two perturbations 
at the critical point could result in the phase diagram in 
FigUlwhich is consistent with our numerical calculations. 
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